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ADAPTIVE WEIGHTING METHOD FOR ORTHOGONAL FREQUENCY 
DIVISION MULTIPLEXED SOFT SYMBOLS USING CHANNEL STATE 

INFORMATION ESTIMATES 

BACKGROUND OF THE INVENTION 

This invention relates to signal processing, and more particularly, to signal 
processing techniques for use in Digital Audio Broadcasting (DAB) systems. 



superior to existing analog broadcasting formats. Both AM and FM DAB signals can be 
transmitted in a hybrid format where the digitally modulated signal coexists with the 
currently broadcast analog AM or FM signal, or in an ^all-digital format without an analog 



because each DAB signal is simultaneously transmitted within the same spectral mask of an 
existing AM or FM channel allocatioiL EBOC promotes economy of spectrum while enabling 
broadcasters to supply digital quality audio to their present base of listeners. Several IBOC 
DAB approaches have been suggested. One such approach, set forth in U. S. Patent No. 
5,588,022, presents a method for simultaneously broadcasting analog and digital signals in a 
standard AM broadcasting channel. Using this approach, an amplitude-modulated radio 
frequency signal having a first frequency spectrum is broadcast. The amplitude-modulated 
radio frequency signal includes a first carrier modulated by an analog program signal. 
Simultaneously, a plurality of digitally-modulated carrier signals are broadcast within a 
bandwidth which encompasses the first frequency spectrum. Each digitally-modulated carrier 
signal is modulated by a portion of a digital program signal. A first group of the digitally- 
modulated carrier signals lies within the first frequency spectrum and is modulated in 
quadrature with the first carrier signal. Second and tliird groups of the digitally-modulated 
carrier signals lie outside of the first firequency 5pectrum;^d are modidated both in-phase and 



Digital Audio Broadcasting is a medium for providing digital-quality audio. 




new spectral allocations 
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in-quadrature with the first carrier signal. Multiple carriers are employed by means of 
orthogonal frequency division multiplexing (OFDM) to bear the commtmicated information. 

FM IBOC broadcasting systems using have been the subject of several United 
States patents including patents no. 5,465,396; 5,315,583; 5,278,844 and 5,278,826. In 
addition, a commonly assigned pending patent application for a "Method and System for 
Simultaneously Broadcasting and Receiving Digital and Analog Signals, by D. Kumar and B. 
Hunsinger, Serial No. 08/294,140, filed July 1994 discloses an FM IBOC DAB system. 

The signals used in Digital Audio Broadcasting are subject to fading and noise 
(interference). Digital Audio Broadcasting receivers may include Viterbi decoders. 
Conventional implemmtations of soft-decision Viterbi decoders rely on constant signal and 
gaussian noise statistics for (near) optimum decoding. Practically these statistics should be 
nearly constant over the path memory of the Viterbi decoder, or tihe span of the interleaver, 
whichever is greater. An interleaver may be used to yield statistical independence of the 
fading statistics of the soft symbols over the path memory of the Viterbi decoder after 
deinterleaving. The remedy for a flat fading channel causing fluctuating signal levels with 
constant noise is well known. However, there exists a need for a signal processing technique 
that can address independently varying signal and noise levels. 

SUMMARY OF THE INVENTION 
A method for weighting orthogonal frequency division multiplexed soft 
symbols is provided, including the steps of receiving a plurality of sub-carriers modulated by 
digital information, filtering the sub-carriers to produce complex soft decision outputs, 
creating a first sequence of the magnitudes of the complex soft decision ou^uts, determining 
the differences between successive samples in the first sequence, creating a second sequence 
of the differences between successive samples in the first sequence, filtering the first and 
second sequences and compensating for differences in effective group delay of the first and 
second sequences to produce third and fourth sequences, using the third and fourth sequences 
to detemiine a plurality of weights, and applying the plurality of weights to the complex soft 
decision outputs. The invention also encompasses receivers that incorporate the above 
method. 

BRIEF DESCRIPTION OF THE DRAWINGS 
Figure 1 is a diagram showing the steps in the method of this invention; and 
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Figure 2 is a functional block diagram showing the mapping and processing of 
bits through a digital audio broadcasting receiver which employs the method of this 
invention. 

DESCRIPTION OF THE PREFERRED EMBODIMENTS 
The present invention will be described in terms of the adaptive weighting 
procedure for coherently detected quadrature phase shift keyed (QPSK) or binary phase shift 
keyed (BPSK) signals under various fading conditions assuming gaussian noise or 
interference statistics. 

For single carrier, flat fading signal, constant noise or interference, the flat 
fading channel hiultiplies the k^ interleaved symbol by a fading factor a^^ , which is assumed 
to be constant over the symbol diuration. The signal sample ou^ut of a matched filter in the 
receiver can be expressed as: 

>^*=«*-^*+«A (1) 
where jc^^ is the transmitted signal component and n^^ is the received noise. The appropriate 

weighting for that soft decision symbol at the receiver is itself , or more practically the 
receiver's estimate . Therefore the soft symbol Sf, entering a deinterleaver prior to the 
Viterbi decoder is presented as: 

^k^^k^yk (2) 
The symbols immediately surroimding the k* symbol are similarly affected by the relatively 
slowly varying fiiding* It is impossible to adequately estimate this fading factor over a single 
s3anbol because it is impossible to separate its statistics from the noise. However the fading 
factor can be estimated statistically using a weighted average over a set of samples with 
correlated fading constants, specifically the symbols that are temporally close to the k* 
symbol. 

Let the signal term jc^ be complex whose real and imaginary components are 
binary valued at ± 1 . For reasonably large SNR, where the noise (or interference) is 
approximated to be zero mean average white Gaussian noise (A WON), the fading fector 
might be estimated as: 
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However, this estimate is not accurate for small SNR since the rectified noise will place a 
lower limit on a,^ . This limitation would reduce erasure effectiveness when the signal level 
fades to zero. The signal level for coherent detection and Gaussian noise can be estimated 
using an exact expression as a fimction of the second and fourth moments of the magnitude of 

5 (i.e. } and } ). These moments are obtained through averaging or lowpass 

filtering of the squared or fourth power of the magnitude of . Assxuning a constant A WG 
noise floor, the weighting factor for the soft decision is equal to the fading factor: 

This equation has been used to success&Uy estimate the signal level in a mobile satellite 
10 receiver. 

For single carrier, flat fading signal, and flat fading interference, the difference 
between the constant signal and the flat fading signal appears as an interference fading fector 
Cf^. It is assumed that the fading of the signal and interference are approximately 
independent, and the interference is approximated as AWGN. The Gaussian interference 
15 assumption is made in the absence of known precise statistics of the interferes Furthermore, 
this assumption becomes less objectionable when we recognize that the Viterbi decoder tends 
to "Gaussianize" the noise or interference summed over the weighted symbols in its path 
memories. The soft-decision ou^ut of the matched filter can be expressed as: 

20 Dividing both sides by cr^ puts the expression in form familiar to the nonfaded interference 
case: 

^ = -^-x,+«, (6) 

Then the appropriate weight for the soft synibol is identified- The soft symbol entering the 
deinterleaver prior to the Viterbi decoder is presented as: 

25 s,--^'y, (7) 

The fading factor can be estimated as in equation (4). An additional estimate of the noise 
variance is required for this weighting. This estimate can be computed as: 

ai=EM'}-al (8) 
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For multiple sub-carriers (OFDM), selective fading signal, and selective 
fading interference, the method of the present invention can be described as follows. OFDM 
sub-carriers comprise a set of orthogonally spaced (in frequency) sub-carriers, each 
modulated synchronously at the symbol rate. QPSK sub-carrier modulation is asstmied here, 
although other linear modulation techniques such as MPSK or QAM can be accommodated 
with modification. OFDM has been shown to be tolerant of multipath fading since the fading 
over any individual sub-carrier bandwidth is considered to be flat (flat fading requires no sub- 
carrier equalization). In a severe selective fading instance, a portion of the OFDM sub- 
carriers may be lost in a null while the remaining sub-carriers are recoverable. In contrast, 
wideband single carrier modulation can experience the same spectral null where significant 
degradation can occur over all the bits if equalization is not employed. 

Each OFDM sub-carrier can be treated as an independent channel with flat 
fading signal and interference; however, this is not optimum. The correlation of the signal 
and noise fading across sub-carriers nearby in frequ^cy can be exploited to improve 
performance. This can be accomplished through filtering the statistics across both time and 
frequency surrounding any particular sub-carrier. The symbol of the n* sub-carrier is 
weighted as: 

s^.-^-y.. ^ (9) 

where is a complex soft symbol resulting firom QPSK demodulation of the k**" symbol of 
the n* sub-carrier. The only difference between equations (9) and (7) is the exploitation of 
the correlation of adjacent sub-carrier statistics in equation (9). The fading factors can be 
estimated after 2-dimensional smoothing (filtering across time and sub-carriers) of the 
following values: 



Compute the symbol weighting factor for coherently detected QPSK as: 



(10) 



where w^^ is the weight to be used for the soft symbol prior to deinterleaving 
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The adaptive weighting procedure for differentially detected QPSK signals 
under various fading conditions will now be described. These results can be extended to 
BPSK with minor modification. 

The goal of the present invention is to maximize the resulting SNR of the 
5 weighted and combined symbols after differential detection into the Viterbi decoder. This 
minimizes the probability of bit error at the output of the Viterbi decoder. The absolute 
optimization of the weighting of the soft decision statistic is a function of the particular 
Viterbi decoder. Therefore the answer is not necessarily unique. Furthermore, since the 
Viterbi decoder tends to "Gaussiahize** the path metrics over the path memories, then this 
1 0 technique tends to maximize the SNR over, the path memory. 

A sunilar expression for the weighting factor can be derived assuming 
Gaussian noise into a differential QPSK detector resultuig in non-Gaussian statistics at the 
ou^ut. The fading factor can be computed as a fimction of fhe statistics of the output of the 
differential detector where we define the soft decision of the form: 
15 5 = (a + «,)-(a-e^*'+W2) (12) 

where denotes the phase information in^}osed between a pair of adjacent sjonbols in the 
differential encoding, and n are the independent noise samples. The fading factor a of the 
adjacent symbols is assumed to be approximately equal. The signal to noise ratio after 
differential detection is easily computed to be: 



20 



a' 



The ideal weighting factor for the post-differentially detected symbols is therefore: 



The first differential approach described here uses statistical estimates of the second and 
fourth moments of the differentially detected symbols to form the weighting factor. The 
25 second and fourth moments are described by the following previously known relationships. 

Then the fading factor can be estimated as : 

a\^'^j2>EiS,f}-^EiS,\'} (16) 
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and the noise can be estimated as: 



a\^^Ei/S,f}-al (17) 
The estimates of equations (16) and (17) are inserted into equation (14) to obtain the weight 

Simulations were performed using adaptive weighting as described in 
equations (14), (16), and (17). Although long-term estimates without fading yielded good 
results, a compromise must be reached between long filter time constants for accurate 
estimation versus short filter time constants needed to track varying statistics due to fading. 

In the Digital Audio Broadcasting (DAB) simulation, the OFDM symbol rate 
of 689.0625 Hz was chosen with a fading bandwidth of 13 Hz. Then the reciprocal of the 
fading bandwidth is about 53 symbols in this case. A filter time constant of 16 symbols was 
chosen since this time constant must be small compared to the fading time. Unfortunately, 
the statistical estimation errors over this short filter time yielded poor performance results for 
the adaptive weighting con:4>ared to what would be possible with perfect statistical 
estimation. Even reducing the fading bandwidtii down to 3 Hz and increasing the filter time 
constant to 64 samples left a significant loss. 

Equations (16) and (17) reveal that, in effect, quantities raised to the fourth 
power are subtracted to yield smaller numbers. This situation is most pronounced when the 
signal and noise powers are approximately equal, resulting in large estimation errors. 
Simulation results support this observation. Therefore another estimator is sought that does 
not rely upon subtraction fourth order statistics. The desired estimation technique should be 
designed to accommodate a fading bandwidth of up to 13 Hz for maximum vehicle speeds in 
the FM band around 1 00 MHz. 

The optimum soft-symbol weight to be applied before differential detection of 
QPSK can be described as a function of time (k index) and OFDM sub-carrier (n index). 
Similar to equation (14), this weight is: 

a. 



(18) 



where a^^, is the fading coeflBcient of the k**' symbol for the n* sub-carrier, and qj,, is the 
corresponding stsuidard deviation of the noise or interference^ both prior to differential 
detection. Notice that the weight of equation (1 8) is the square root of equation (14). This is 
a result of the reasonable assumption that the weight changes slowly over the symbol-pair 
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time used in the differential detection. In effect, the differential detection squares the 
predetection weight of equation (18), which would result in equation (14). A method for 
improving the statistical estimates of equation (18) is sought 

Practical methods for estimating CSI and weights using pre~differentially 
detected soft-symbols and weight also applied to the soft decision symbol prior to differocitial 
detection will now be discussed. 

For moderate to high SNR, the weight of equation (18) can be conveniently 
approximated by: 

= / , (19) 

where simple statistical measurements were used to estimate o-^. However, simulation 
confirmed that this weight estimate performed poorly during times when the SNR was very 
low due to fading interference. For example, the optimum weight would have suppressed the 
noisy samples more than the high SNR approximation to the weight. Therefore, another 
approximation was sought which would estimate CSI statistics over a large SNR range. 
Furthermore the estimate should not be sensitive to a gaussian noise or interference 
assumption, and should be estimated with sufficient accuracy in a time (filter time constant) 
significantly less than the reciprocal of the fading bandwidth. 

A simple and robust estimation technique evolved after simulation and some 
experimentation. This estimation technique approximates the previously-defined weight 
expressions, but uses lower-order statistical approximations. This technique is described in 
the following 4 steps. 

1. Create a sequence Vj^^ for eadi QPSK sub-carrier consisting of the magnitudes of the 
complex soft decision outputs s^^^, from the matched filter for the n* sub-carrier. 

(20) 

2. Create a sequence df^^^ consisting of the differences of successive time samples of 

3. Filter the sequences v^^^ and c/^^, using second-order digital IIR filters, then 
compensate for any differences in effective group delay to yield sequences ftltVi,^ and 



wo 00/16528 



- 9 



PCT/US99/20275 



filtdf^^. The time constant for the filtv^^ filter should be somewhat smaller than the 
reciprocal of the fading bandwidth, while the time constant for the filtdf^j, filter can be 
somewhat larger. These sequences are representative (approximately proportional) of the 
local mean and standard deviation of the sequence v^^^, . 

4. The sequence of weights for the soft decisions for each sub-carrier to be applied prior 
to differential detection is defined as 

1 

(22) 

To prevent numerical overflow, check to ensure that filt^k^ >1^- fil^dg^^ in equation (22); 
otherwise, set tibe weight to zero. Simulation results verified that this weight yields good 
performance imder a variety of channel impairments with fading and interference^ 

The values of ftltd^^ and >!/<Vjt^ are estimated using filtering techniques 
described next. Filtering is performed first for each sub-canier at the k"" symbol instant in 
time. Then the rows of flltdf^^^ and filtVf^^ are simply updated across the N sub-carriers. 
Equation (23) filters the sequences v^^^, with a time dejay of approximately 16 symbols, and 
equation (24) filters the sequences d^^ with a time delay of approximately 64 symbols. Both 
filters have a zero firequency gain of nearly unity. 

su}^ = !^ (23) 

512 

subd.^^ ^ ^ ^ (24) 

*,« 8192 

Additional filtering is performed across the N sub-carriers. Smoothing the 

estimates across the N sub-carriers requires 3 passes of a simple nR filter. The first pass sets 

the appropriate initial condition of the filter, but does not update the estimates. The direction 

of the second pass is reversed firom the first, while the itoird pass is reversed again. This 

results in an approximately symmetric (linear phase) filter characteristic which is desirable 

for providing the estimates on the center carrier. Although it is impossible to provide this 

symmetric filtering for the sub-carriers at each end of the band, the impulse response *^ails" 

are folded back into the active sub-carriers. 
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The first pass across the sub-carriers sets the initial values of filtv^^^^ and 
filtd^,^ without replacing the time-filtered values for each sub-canrier. The time index k is 
ignored here smce it is understood that the filtering over the sub-carriers is performed over 
each OFDM symbol. 

filtv^,, c=(l->ff)-yi/A;^., -^fi subv, ; 

flltd^,, <=i{\-fi) filtd^^, ^fi^subd, (25) 

The second pass smoothes the values across the filtered estimates for each sub- 
carrier, subv and subd . 

filtd, <;=(l^JJ) flltd^, ^P subd,\ (26) 

/i = iV--2,A^-3,...0 

The third pass smoothes flie firequency values again to achieve a nearly 
sjrmmetrical impulse response (except for the sub-carriers near the endpoints). 

filtd, <^(\-m^filtd^x ^P'filtdnl (27) 

The resulting filtered values for ftltv and filtd computed in equations (26) and (27) are used 
in equation (22) at each OFDM symbol time to yield the appropriate weight for each soft 
symbol prior to differential detection, but after matched filtering, in the receiver. 

Figure 1 is a digram showing the steps of the method of this invention. Block 
10 shows the step of receiving a plurality of sub-carriers modulated by digital information. 
The sub-carriers aie filtered as shown in block 12 to produce complex soft decision outputs, 
these outputs are used in block 14 to create a first sequence of the magnitudes of the complex 
soft decision outputs. Block 16 shows that the differences between successive samples in the 
first sequence are determined. A second sequence of the differences between successive 
samples in the first sequence is created as shown in block 18. The first and second sequences 
are filtered as in block 20 and compensated for differences in effective group delay of the first 
and second sequences as in block 22 to produce third and fourth sequences in block 24. The 
third and fourth sequences are used to determine a plurality of weights as shown in block 26, 
and tlie plurality of weights are applied to the complex soft, decision outputs as shown in 
block 28. 
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The above discussion relates to dififerentially detected QPSK. The derivation 
for the weight using coherent detection of QPSK or (BPSK) would be similar to the technique 
discussed above for differentially detected QPSK. The only modification should be in the 
expression presented in equation (22)- It can be predicted that the weight should be of the 
form : 



where c is a constant to be defined by empirical methods. 

Techniques for estimating tiie optimum soft-decision weight for QPSK 
symbols prior to Viterbi decoding were described. These techniques apply to coherent and 
differential detection of single or multi-cairier (OFDM) QPSK, with and wi&out multipath 
fading of the signal of interest or the interferes The fading cases can necessitate a 
compromise between accurate CSI estimation and agility of the CSI to track the fading signal 
or noise components. 

The application of soft-decision weighting for an OFDM system with 
independently faded signal and noise (interference) can improve BER performance. This 
weight is applied to the soft symbols prior to deinterleaving and Viterbi decoding. The 
optimal weight for each soft symbol over time and sub-carriers is estimated through an 
expression using filtered statistical estimates of channel state information (CSI) about the 
signal and noise components of the received symbols. The above CSI estimation and 
weighting techniques can apply to both coherent and differentially detected symbols, and to 
single and multi-carrier (OFDM) modulation in the presence of multipath fading and colored 
noise statistics. 

Figure 2 is a fimctional block diagram showing the mapping and processing of 
bits through a of a portion of an FM receiver that operates in accordance with the method of 
the invention. A plurality of OFDM carriers 30 are received and converted to bit streams on 
lines 32 by receiver circuit 34. Circuit 34 includes a digitizer, carrier synchronization, 
symbol synchronization, and matched filters all operating in accordance with well known 
techniques to produce the bit streams on lines 36. Block 38 represents the channel state 
estimates and weighting processes performed according to Figxire 1 to produce weighted bit 
streams on lines 40. Block 42 shows that the bit streams are deallocated from the carriers and 
delivered to a deinterleaver 44. The output of the deinterleaver is multiplexed to a single bit 
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stream and passed to a Viterbi decoder 46 decodes the single bit stream. Soft-decision 
Viterbi decoding with (near) optimum soft-decision weighting for maximum ratio combining 
(MRC) for differentially detected QPSK sub-carrier symbols is employed to minimize losses 
over the channel- The output of the Viterbi decoder is subject to additional signal processing, 
which is not part of the present invention, and passed to output stages, as illustrated by block 
48, to produce Ifae desired output fi-om the receiver. 

While the present mvention has been described in terms of what are at present 
believed to be its preferred embodiments, it will be apparent to those skilled m the art that 
various changes miay be made to the embodiments described above without d^arting fix>m 
flie scope of the invention as set forth in the foUowmg claims. 
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What is claimed is: 

1. A method for weighting orthogonal frequency division multiplexed 
soft symbols, said method comprising the steps of: 

receiving a plurality of sub-carriers modulated by digital information; 

filtering tiie sub-cairiers to produce complex soft decision oulputs; 

creating a first sequence of the magnitudes of said complex soft decision 

outputs; 

determining flie differences between successive samples in said first sequence; 
creating a second sequence of the differences between successive samples in 
said first sequence; 

filtering said first and second sequences to produce third and fourth 

sequences; 

using said third and foxirth sequences to determine a plurality of weights; and 
applying said plurality of weights to said complex soft decision outputs. 

2. The method of claim 1, wherein said sub-carriers are modulated using 
quadrature phase shift keymg. 

3. The method of claim 1, wherein said second and third sequences are 
filtered using a second order digital infinite impulse response filter. 

4. The method of claim 1 , wherein: 

said third sequence has a time constant smaller than the reciprocal of the 
fading bandwidth of plurality of sub-carriers; and 

said fourth sequence has a time constant larger than the reciprocal of the 
fading bandwidth of plurality of sub-carriers. 

5. The method of claim l,wherem the step of filtering said first and 
second sequences fiirther comprises the step of compensating for differences in effective 
group delay of said first and second sequences to produce third and fourth sequences. 
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6. The method of claim 1, wherein said plurality of weights are 
detenxiined using the formula: 

1 



w.. = 7 : TT\ 



wherein Wk,n represents said plurality of weights, represents said third 

5 sequence, and flltdk^n represents said fourth sequence. 

7. The method of claim 6, A^eiein i£filtvk,n > L5 ^filtdk^n. then is set 

to zero. 

8. A radio receiver comprising: 

means for receiving a plurality of sub-caniers modulated by digital 

10 information; 

means for filtering the sub-caniers to produce complex soft decision outputs; 
means for creatmg a first sequence of the magnitudes of said complex soft 
decision outputs; 

means for determining the differences between successive samples in said first 

15 sequence; 

means for creatmg a second sequence of the differences between successive 

samples in said first sequence; 

means for filtering said first and second sequences to produce third and fourth 

sequences; 

20 means for using said thurd and fourfli sequences to determine a plurality of 

weights; and 

means for applying said plurality of weights to said complex soft decision 

outputs. 

9. The receiver of claim 8, wherein said sub-carriers are modulated using 

25 quadrature phase shift keying. 

10. The receiver of claim 8, wherein said second and third sequences are 
filtered using a second order digital infinite impulse response filter. 

11. The receiver of cl2iim 8, wherein: 

said third sequence has a time constant smaller than the reciprocal of the 
30 fading bandwidth ofplurality of sub-carriers; and 
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said fourth sequence has a time constant larger than the reciprocal of the 
fading bandwidth of plurality of sub-carriers. 

12. The receiver of claim 8, wherein the means for filtering said first and 
second sequences further comprises means for compensating for differences in effective 
group delay of said first and second sequences to produce third and fourth sequences. 

13. The receiver of claim 8, wherein said plurality of weights are 
determined using the formula: 

1 



wherein wjt,„ represents said plurality of weights, ^//Vit,rt represents said third 
10 sequence, and filtdu,„ represents said fourth sequence. 

14. The receiver of claim 13, wherein if Jzltv^n > I S •filtdk,n> then wjtn is 

set to zero. 
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